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NOTATION 


Total beam cross sectional area. 
Beam Flange Area. 


Cross sectional area of beam web 
and top flanges. 


GGerEtoTent Macrix Of COnstraint 
equation. 


Distance between centroids. 

Constants of constraint function. 

Vector of the basis. 

Right hand side of constraint relation. 
Constant of constraint function, 

Length of beam overlap. 

Coefficient matrix of objective function. 
Bracket leg length. 

Objective function, 


Lengths of rigid sections due to beam 
overlap. 


Beam element weight per unit length. 
Beam depth. 

Beam element length. 

Plastic limit moment. 

Vector associated with extreme point. 
Overlapping web depths. 

Cross-tie height above base line. 


Side tank width. 
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NOTATION (cont. ) 


Rigid extension due to bracket 
influence. 


Normal velocity. 
In line velocity. 


Variable of linear programming 
problem. 


Iteration solution value of objective 
FUNG E LON. 


Plastic section modulus. 
Multiplier vector. 
Material density. 


Material yield stress. 





1. INTRODUCTION 


1.1 Preface 

The design problem of transverse web frame 
structures in ships is receiving an increasing amount of 
interest with the advent of mammoth tankers and bulk 
Carriers. Transverse framing accounts for a sizable 
portion of the total structural weight of such ships. A 
substantial weight savings can result from the proper 
design and optimization of the repeating transverse 
frames in the prismatic mid section of these large vessels. 
Great effort is usually made in the design phase to make 
these structures as efficient as possible. Modern day 
design efforts have been concentrated in two methods of 
elastic design. 

The finite element method is the most popular of 
the modern elastic ship design methods. A number of 
review papers [1 - 4] have discussed the application of 
this technique to the design of ship structures in general. 
Smith & Woodhead [1] discussed specific application of the 
method to web frames. 

The second elastic design method uses beam theory 
and was first discussed in depth by Lund [5]. Others have 
also written on the topic [4 - 14]. 

Both of the elastic design methods involve 


systems of nonlinear equations. Design optimization can 





become complex and expensive in terms of computation time. 
The nonlinearities require the use of optimization 
techniques such as nonlinear programming, Powell's direct 
search method or the gradient method. 

The elastic design approach and the nonlinear 
Optimization methods combine to form a laborous design 
task with a trusted but possibly expensive and inefficient 
product. 

The purpose of this study is to examine the use 
of another design method, namely plastic limit design, 
in the design of transverse web frames in ships. 

1.2 Historical Review 

A search of existing literature reveals a wealth 
of information on the use of plasticity design methods and 
specifically the use of the plastic limit design theorems 
dating back several decades [15 - 22]. The limit design 
methods are common place and vital in nearly all sectors 
of structural design except naval architecture and ocean 
engineering. In papers on marine structures where the use 
of limit design is suggested [23 - 25], it is often from 
the civil engineering point of view with examples taken 
from standard civil engineering texts. 

Jones [24, 25] uses the limit theorems in the 
analysis of tanker web frames in two separate papers. To 
date, his are the only papers known to discuss in detail 


the use of the limit theorems on marine structures loaded 





by various conditions of hydrostatic loading. 





2. METHOD OF CALCULATION 


Pee ie Sling aby Design Upper Bound Theorem 


In structures composed of an elastic, perfectly 
plastic material, one of the primary functions of the 
structure is that it should carry a given set of loads 
while in equilibrium without flow occurring in the 
structure. The structural design problem is referred to 
as the limit design problem if this is the only function 
of the structure [22]. An alternative formulation of the 
limit design problem states that if under all possible 
failure modes the internal energy dissipated during 
kinematically admissible, small virtual displacements 
exceeds the external work rate associated with these 
displacements then the structure is safe [7]. This 
latter statement can be used to formulate an upper bound 
on the system plastic moments as a programming problem. 

The limit design upper bound theorem is the 
primary tool to this analysis. The desired output of the 
limit design problem is the required weight (i.e. cost) of 
individual structural elements and the total system. 

The output of the limit design process is the 
plastic moment, therefore, it is necessary to establish 
some relationship between weight per unit length g and 
plastic moment M,. Horne [21] states that for materials 


of given yield stress and geometrically similar section, 
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g «= M273, For Universal Beam Sections the relationship 
has been shown to closely approximate g « M*®. More 
generally g « M™", Horne also states that the relative 
plastic moments for absolute minimum weight in a given 
structure are only very slightly affected by the exact 
value chosen for the index "n". It will be shown that 
solutions are greatly simplified by maintaining linear 
relationships between variables. It, therefore, is 
assumed that the weight per unit length of an element 
is given by 

g=AL+BLM, 
where A and B are constants and L is the element length. 
The total system weight of the structure is then given 
by 

Glee eet) Be bMS 

The methods of plastic limit design cannot 

utilize the principle of superposition of stresses and 
deflections due to individual loads, as is possible in 
design where the structure remains in the Sieecie range 
[15]. A paradox of the limit theorem, however, is that 
when a collapse mode of a structure involves only part 
of the structure, the collapse load of the plastic limit 
moment can be expressed independently of the characteristics 
of the remainder of the structure [18]. In other words, 
in formulating the limit design problem, only those loads 


associated with a given collapse mechanism need be 
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considered with that mechanism. This premise allows the 
energy equations associated with each mechanism to be 
constructed independently. 

Further assumptions and idealizations are 
required to complete the formulation of the limit 
design problem. Material behavior is idealized to be 
elastic, perfectly plastic (i.e. free of strain hardening). 
Strains are assumed to be small. Jones has discussed 
the importance of these assumptions [26]. 

Finally, it is assumed that all frames are 
plane frameworks of beams, in which transverse shear and 
axial forces do not influence the plastic yielding of the 
beam. Symonds and Neal [15] stated that shear forces 
can be neglected except in cases where the ratio of span 
length to beam depth is less than approximately four to 
one. Also, the effect of an axial force on the fully 
plastic moment can be neglected when the force is less 
than about one fourth of the compressive yield force 
of the member. 

To illustrate the general proceedures of 
formulating a limit design problem, the simple civil 
engineering portal frame of Figure 1, will be examined. 
Application to the more sophisticated ship frames with 
hydrostatic loads is discusses in Reference [25] in the 
appendices and in sections 2.2 to 2.4. The frame is 


restrained vertically as well as horizontally at its 
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base and the supports are clamped against rotation. 

The frame is externally determinate and internally 
indeterminate with r = 3 redundancies. A discussion of 
redundancies may be found in [27]. 

The frame has n = 8 critical sections, or 
locations where plastic hinges can form. Plastic hinges 
may form at the ends of a span, directly under a 
concentrated load or within the span of a member 
subjected to a distributed load. 

Only r of the n critical bending moments can 
be defined independently of the equilibrium equations 
[24]. Therefore, there exist n - r or 5 independent 
relations between the bending moments. These relationships 
are the system equilibrium equations and can be obtained 
by utilizing the principle of virtual work. Each 
equilibrium equation is associated with a kinematically 
admissible fundamental mechanism. Fundamental or 
elementary mechanisms are classified by Symonds and Neal 
as beam, panel, sway and joint rotation. Joint 
rotations are also classified fictitious. The 5 
mechanisms of the sample problem are illustrated in 
Figure 2. For all 5 failure mechanisms to safely 


carry the prescribed loads we must have 


M, - 2Mg + Mg 2 2 (2-1) 
M,; - 4M) + 3M, > 6 (2.2) 
M, - M3 + Mo - Mg 2 6 (253) 
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WS ily ele (2.4) 


Me - M5 > 0 (255) 


7, 
where values of 6, angular velocity have been divided 

from both sides of the inequalities. Minus signs denote 
hinge rotations in which the internal fibers of the 

hinges are in tension. This convention is merely a 

method of keeping track of the direction of rotation 

and is not intended to imply negative internal energy 
dissipation. 

The fundamental sway mechanism was eliminated 
from consideration in this study. This was not because 
of any assumptions that sway would not cause collapse 
in ship web frames, but rather that sway could not be 
induced by any of the symmetrical loading models used 
by the elastic designs which are used as comparison 
baselines for this study. 

The fundamental mechanisms are not the only 
mechanisms which must be examined. It is possible for 
linear combinations of the fundamental meonenitene to 
control the design [28]. This is true where the combination 
of mechanisms allows a reduction or elimination of 
certain plastic hinges reducing the total internal 
energy dissipation while the combined external work 
rate remains unchanged. It is here that the signs 
of the various terms in the equilibrium equations become 


helpful bookkeeping aids. In combining the fundamental 


aA 


mechanisms of the example the fictitious joint rotations 


are used to derive 


M, = 4M, + 3M, > 6 (2a) 
M) - My + Mo - Mg > 6 (2.8) 


also equations (2.2), (2.3) and (2.5) can be combined to 
give 
4M, - 4M, + 3M7 - 3Mg 2 24 (229) 
or an auat hol = oho 224 (210) 
Assume, now, that the structure will be 
developed with two separate beam sizes. The two vertical 
elements will be identical in cross section with plastic 
moment M_ and the horizontal element is allowed to differ 


from the verticals with plastic moment M Now equations 


De 
(22) 2 thirougis (2.s)mand, (276) through: (2.10): can be 
written in terms of two variables with the sign 


designators removed as follows 


4M, 2 2 (2.1la) 
8M. 2 6 : (22224) 
4M, > 6 (2.3a) 
2M,+ 2M, > 2 (2.6a) 
Soe Me 6 (2. 7a) 
2M,+2M, 2 6 (2.8a) 
14M, > 24 (2.9a) 
11 M, + 3 My 2 24 (2.10a) 


Note that the fictitious joint rotations can not physically 


SS 


occur alone and therefore appear in this final formulation 
only in the formation of the combined mechanisms. The 
minimum weight design problem is then to minimize 


Minf=- 2LM +L, M 
a a 


| oad ©, 


or Min f= 8M, +4™M (2510) 


b 
subject to expressions (2.la) to (2.3a) and (2.6a) to 
(2.10a), where the function (2.11) is the design minimum 
weight function. 

It is here that the assumption of a linear 
weight and moment relationship discussed on page 4 
becomes important. In the present linear form the 
problem of the minimum weight design becomes a relatively 
Simple linear programming problem in two variables, Me 
and Mie where the object function is (2.11) and expressions 
(2.la) to (2.3a) and (2.6a) to (2.10a) are the constraining 
functions. 

This linear programming problem will be solved 
for M. and M, in section 3.1. 
2.2 Rigid Bracket Region Determination 

In the limit design problem of very large and 
deep ships, with deep web frame structures, the effect 
of web overlap and bracket influence becomes significant. 
It is customary [5, 14] in elastic design of ship frame 
structures using beam theory to assume that a certain 


portion of brackets and overlapping webs constitute a 


rigid region on the ends of each span as shown in Figure 
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3. This method called the "Span Point Method" has been 
shown by finite element calculations to slightly over 
extimate the extent of the rigid regions. 
To calculate the length of the rigid section 
in a given element Lund [6] uses the expressions 
fr = Srr f/f sina 
for = Sy / Sina 


where fs and for are the lengths of a rigid section in 


beams I and II respectively and S_ and §S are the 


I Il 
overlapping web depths associated with each beam. The 
value a is the angle of intersection between the two 
beams aS shown in Figure 4a. The angle is measured at 
the vertex of the overlapping webs on the beam ends and 
is formed by the lines parallel to the beam axes through 
that vertex. It is evident that in the most common case 
of normal beams butting flush against each other, the 
rigid length influence on each beam of the intersection 
is the depth of the adjacent beam. Where three or more 
beams tie in at the same point, the beams mea the greatest 
depth govern in this effect. In this case the bracket 
influence has been neglected. 


To compute the influence of each bracket leg, 


the approximation 


b 
i 


is used [6]. As is shown in figure 4b., ay is the 


t, = 4d, /{(h/d,) + 1} 


computed rigid length in the bracket, qj, is the bracket 
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leg length and h is the beam depth. The total rigid 
length of one end of a given beam is the sum of the 
effects of the bracket and the beam overlap Cs, 
Figure 4b. 

In this study, when bracket influence is 
considered, intersecting beams will be normal or nearly 
normal and it is assumed that c, = f,. 

It should be noted here that the preceeding 
expressions are based on known values of geometry which 
are not known at the start of the design process. 
Iteration from an assumed starting value is reguired. 

The rigid region assumption takes on special 
Significance in the limit design problem. The effect is 
to indirectly specify regions within the structural system 
in which plastic hinges may not form. This forces span 
ends toward their midpoints while the total external work 
rate part of the virtual work relation may or may not 
remain unchanged, depending upon the nature of the 
particular mechanism under consideration. In any event, 
the formulation of the virtual work relation becomes 
Significantly more complex in all but the simplest beam 
mechanism, especially in the external work part of the 
expression. The variability of rigid length values 
through a successive iteration steps compounds the 
problem. Lund [5] side steps this problem by calculating 


the rigid extensions only once at the problem beginning 
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and holding them constant thereafter. 

Rigid section examples examined in this study 
were Calculated by hand with rigid extensions assumed to 
be the same as used by Lund. 


2.3 Combined Mechanisms and Kinematic Admissibility 





Two final considerations are very important 
in the plastic limit design calculations performed in 
this study. These are limitations which can restrict 
combining mechanisms as was illustrated in section 2.1 
and conflict with the requirement that all mechanisms 
be kinematically admissible, 

The upper bound theorem of limit design requires 
that all velocity fields be kinematically admissible. 
This is a requirement that deserves special attention in 
formulating the limit design problem. For the velocity 
field in a plastic collapse mechanism to be kinematically 
admissible the velocity in the direction normal to the 
axes of the elements, u must be much larger than in line 
velocity, w. Figure 5 illustrates that ale slight 
increases of the beam lengths are required to move point 
b a substantial distance laterally through beam rotation 
to point b' in the pictured mechanism. 

Martin [22] states that the hinges of joint 
rotation mechanisms must form at the joint and not some 
finite distance out. This point is generally overlooked 


in the literature, but is of utmost importance if the 
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joint mechanism is to be utilized in forming combined 
mechanisms. The source of this problem is the requirement 
that the displacement fields of all limit design 
mechanisms be kinematically admissible. If a multi- 
element joint is rotated and if its hinges are located 
anywhere but at the center of rotation, kinematic 
admissibility as defined previously will be violated. 
Plastic hinges in addition to the joint mechanism 
hinges will be required before the mechanism can rotate. 
The assumption of a rigid bracket region 
conflicts with the use of joint rotation mechanisms 
Since the rigid regions permit hinges at their extremes 
away from the rotation center. Therefore, while the 
rigid joints can certainly be used as components of 
mechanisms, they cannot be considered as independent 
mechanisms. This all but eliminates linear combinations 
of mechanisms as discussed in section 2.1, if joint 
rotations are required to form the combination. When 
combined mechanisms are formed they require the 
formation of a very large number of plastic hinges 
with correspondingly high values of internal energy 
dissipation, reducing the probability that such a 
mechanism would contribute to collapse. In this study, 
combined mechanisms were considered only when they 
consisted of seven or less plastic hinges. 


The use of the rigid bracket region assumption 
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can be one final source of trouble. If the rigid 
extensions are large and angles between hinge points are 
large, the proposed mechanism may be kinematically 
inadmissible. In Figure 6 which approximates the rigid 
bracket in a bilge region, u and w are of the same 

order of magnitude and the mechanism is inadmissible 

as 1s. Were hinge c moved farther to the right the 
mechanism could become admissible. How far the hinge 
must be moved before an admissible mechanism will be 
formed is a question of judgement. For this study, it 
was generally required that the height of the elevated 
hinge, h, be less than twice the distance from the 
corner below it to the next plastic hinge for the 
mechanism to be considered kinematically admissible. 

It 1s assumed that if this requirement is met, the value 
of a will be relatively small. 

The rigid extension assumption is likely to be 
of greatest importance in tanker web frames which do not 
have cross ties. These frames will have — deep webs 
and correspondingly large overlap regions. One frame 
of this type was examined. In frames which have cross 
ties, the web depths will be less. In the transverse 
frames which had cross ties, the rigid extension 
assumption was not used. In frames of this type, plastic 


hinges were permitted anywhere in the beam element. 
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2.4 Summary 
In this study, the upper bound theory of limit 


design is used to develop a system of linear relations 
which must be satisfied if the design structure is to 
safely carry the prescribed loads. A linear weighting 
or object function is formulated based on the lengths of 
the respective structural elements. These linear 
relations combine to form a linear programming problem 
which may be solved to find the minimum required plastic 
moments of the various structural elements in the system. 

It is assumed that the plane frameworks of 
beams are made of materials with elastic, perfectly 
plastic moment curvature relations. It is also assumed 
that transverse shear and axial forces do not influence 
the plastic yielding of the beams. 

Beam, panel and fictitious joint rotation 
mechanisms are considered in the formulation of linear, 
principal of virtual work, constraint relations as well 
as many combined mechanisms. Combined mechaniiane which 
are very complex and fundamental sway mechanisms are 
not considered. 

The assumption of rigid extensions in the 
area of beam overlap and brackets used in some elastic 
designs is considered in one case. This assumption 
disqualifies the use of fictious joint mechanisms in 


combining fundamental mechanisms on grounds of kinematic 
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inadmissibility. Mechanisms involving these rigid 
regions must be carefully checked for kinematic 
admissibility. The rigid extensions, where used, are 


calculated only once at the start of the design 


sequence. 
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36,” OF LIMEZA TION TECHNIOUE 


3.1 Linear Programming 
It was stated in the preceeding chapter that 
the linear equilibrium constraint relations and the 
minimum weight object formed a linear programming 
problem. 
The general linear programming problem can 
be described as follows 
Given a set of i linear inequalities or 
equations in j variables, we wish to find 
non-negative values of these variables 
which will satisfy the constraints and 
maximize or minimize some linear function 
of the variables [30]. 


Mathematically this can be written minimize 


(or maximize) the linear objective function, 


} 
f = ean or (3615) 
j=1 JJ 
subject to the linear equality constraints 

) a. ~R, = Dp. i prove sue (3) 

j=1 1) J Bk 
and 
eas? 2 

x2 ( ) 

The values of a:., b., c. are assumed to be known 


a2 baie Gamat 
constants. Other notations, expecially vector forms, 


also appear in the literature. The use of inequalities 


is not entirely prohibited in equation (3.2) [29]. 
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the use of slack variables permits the inequality 
N 
eee D (3.4) 
J 


to be expressed as 


N 
) a4 5%; eon (3.5) 


The use of a plus sign in front of the slack 
variable would signify that the relation was of the 
less than (<) form. 

Any set of he which satisfies the constraints 
(3.2) is called a solution to the linear programming 
problem. Any solution which satisfies the non-negativity 
restrictions of equation (3.3) is called a feasible 
solution. Any feasible solution which optimizes the 
objective function is called an optimal feasible solution. 
Generally, there will be an infinite number of feasible 
solutions to a linear programming problem. From all 
of these feasible solutions, one must be found which 
optimizes the objective function. 

Linear programming problems which involve 
only two variables can be solved graphically. The 
graphic solution to these two variable problems can 
provide a great deal of insight into the more general 
linear programming case with any number of variables. 


In Chapter 2, a simple portal frame was 
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examined and a system of linear constraining inequalities 
was developed which gave upper bounds on the plastic 
limit moments required to make the structure safe. 
These inequalities in two variables with the minimum 
weight objective function form a linear programming 
problem which is illustrated graphically in Figure 7. 

The solid lines represent the various constraint 
functions. The dotted lines are the parallel lines of 
constant weight of the object function. The solid 
lines with cross-hatching attached bound the region of 
feasible solutions (or the acceptable design region). 
This is the boundary of the region where all inequality 
constraint relations are satisfied. The optimum design 
is found at the point, or points, where the object 
function 

Min f = 8M, pecdule 
has its lowest value while just touching the boundary of 
feasible solutions. The figure shows that the variables 
of the optimal solution are approximately M. = 1.8 and 
M, = 1.1. To find the exact values, simultaneous 
solution of the equations of the two intersecting lines, 
equations (2.8a) and (2.10a) is required. Restated, 
these are 

I, M, + 2 Mp = 6 

11 MA + 3M = 24 


b 


Solution of these equations gives exact values 


29.6 = 





of M = 1,675 ang M —= 1.125. °These values represent 


b 
exact, upper bound solutions to the minimum values of 
plastic limit moment required to make the design 
structure of Figure 1, safe under the prescribed loads. 
Note that no factors of safety have been applied at this 
point. 

In this case the optimum solution was defined 
to be at single vertex formed by two boundary lines. 
It is also possible that the objective function could 
have passed through two vertices as well as all points 
on the line between them. In cases of this type, it 
is said that the linear programming problem has 
alternative optima. This condition occurs when the 
objective one of the extreme points is optimal. 
Unfortunately these methods do not tell which one Pear 

Other methods must be used to solve linear 
programming problems. There is presently no method 
available which will find an optimal solution to a 
linear programming problem in a single step. All of 
the proceedures for solving these problems are, therefore, 
iterative. If the number of constraints is large, it 
becomes very difficult to find any feasible solution 
and still more difficult to find an optimal solution. 
3.2 The Revised Simplex Method 

The revised simplex method of solving linear 


programming problems is a revision of the simplex method 
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derived by Dantzig in 1948. The term "simplex" has 
nothing to do with the technigue as it is now used. 
To properly introduce the revised simplex method it 
is helpful to start with a short discussion of the 
simplex method. Both technigues are algebraic 
iterative procedures which will solve exactly any 
linear programming problem in a finite number of steps, 
or give an indication that the solution is unbounded. 
There are a finite number of extreme points and if 
there is an optimal solution, one of the extreme 
points is optimal. The simplex method is a technique 
which moves from a given extreme point to another 
adjacent and more optimal extreme point, along the 
boundary of the feasible solution region. Of all 
adjacent extreme points, the one chosen is the one 
that gives the greatest decrease (or increase) in the 
object function. At each iteration, the simplex 
method defines whether or not the associated extreme 
point is optimal and if not, it defines the next 
extreme point to be examined. 

Complete discussions on the simplex and revised 
Simplex methods of linear program are by nature quite 
lengthy. It is not the purpose of this study to 
examine in depth the mathematical rational for either 
of these methods. Detailed information on theory, 


development and application can be found in References 
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The simplex method has two phases. The first 
phase finds a first basic feasible solution which 
satisfies the condition that the m constraints intersect 
and that there are only m nonzero variables ae These 
“Ea variables are called the variables of the basis. 

The remaining (N - m) of the variables a are equal 

to zero and called non basic. At the end of the first 
phase the results are structured in standard canonical 
form. In order to satisfy inequality (3.3) all values 
of b. must be greater than zero [8]. In practice this 
may require the introduction of artificial variables. 

When the canonical form of the first phase is 
obtained, the second phase begins. Given any basic 
feasible solution an iteration step 1s made toward an 
optimal feasible solution by changing a single vector 
in the basis. In the revised simplex method, the 
criteria used to determine the vector to enter and leave 
the basis are identical with the simplex method. Changes 
in the basic feasible solution always lead back to the 
canonical form with the variables x renumbered. The 
essential element that allows progression from each 
basic feasible solution is the explicit knowledge of the 
representation of the vectors not in the current basis in 
terms of the basis vectors. 


The main difference between the simplex method 
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and the revised procedure is that in the simplex 

method all elements of the simplex tableau are transformed 
by means of elimination formulas, while in the revised 
simplex method only the elements of an inverse matrix 

need by transformed using these formulas. 

Since its development the revised simplex 
method and in particular a variation which employs the 
product form on the inverse, has been used in the larger 
high-speed computers [29]. 

There are several reasons for the current 
widespread use of the simplex method. For problems 
with sparse coefficient matrices, computation is reduced. 
The revised procedure always deals with original 
coefficients. The amount of new information required 
to be stored is less than in the simplex method, since 
only the inverse and the solution vector must be stored 
and not the complete simplex tableau as in the simplex 
method. 

To formulate the computational procedure for 
the revised simplex method equations (3.1) through (3.3) 
will first be restated in vector form. In the general 
linear programming problem a minimum value of the object 
function 

cx (3275) 
is sought subject to 


AX=b (3.6) 
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and 
X > 0 (a) 
Fe~ting the matrix B of the basis of m-dimenslonal vectors 


Gorrespond to the first m vectors of A, such that 


B Xo = J) (323) 
X,) 2 0 (3.9) 
where Xo = (X49: Xo Qr sees ae) is the solution vector. 


A basic feasible solution is given by 


Xp = Bt b (3.10) 


The linear combinations of all the vectors of 
Piiacermc Of Becan pe Found by 

x. = BTL p. pyle Daten (2) 
where x. is a column vector and P3 1s a set of linearly 
independent vectors associated with each extreme point 
and can be expressed 

P, = X14 Pi + X54 P. eo ee Xm Pi. (3 gaez } 

The initial basic feasible solution is found 
by setting X 5 = a 

For a given iteration solution the value of 


the objective function is 


¥ == + * + @@e@@ @ + * 
Z Cy aie Cy x C_ ae 
ata, seep el, toe} 


where Cc; are cost coefficients of those vectors in the 


basis. Converting to vector form 2. = Cy X. where Co 


2 
= (cy, Co sees cy) is a row vector. 


From equation (3.11) a can be written in the 


ee a 


foOmm 


A (3254) 


J 
where T= fp 


Fae) 


J 
-1 mis) 
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nm is an m-dimensional row vector called the multiplier 
vector. The individual values Ts in the vector are 
called the simplex multipliers. Given the value of 

m for a given basis B, the corresponding ai can be 
calculated. A vector P 5 not in the basis is evaluated 
by computing 7 Ps - Cc, = Ze — C5. 

The information needed to proceed from one 
feasible solution to the next improved feasible solution 
Poetne Original data consisting Of A, b and ¢ and 
equations (3.10), (3.11) and (3.15), 1f£ for each feasible 
basis B, Bl is known explicitly [29]. This concept is 
fundamental to the revised simplex method of solving 
linear programming problems. 

The revised simplex linear programming procedure 
has been applied by large numbers of users to a wide 
variety of problems with efficient results. Most often, 
it can be treated as a "black box" operation with only 
general knowledge on the part of the user. It forms a 
part of most computer application libraries and is 
available at little or no increase in cost over the basic 
computer usage fees [8]. 

The subroutine called in this study was the 


Mathematical Programming System Extended (M.P.S.X.) as 


ee 


installed on the IBM 370 Model 168 computer. Specifics 
on the use of the subroutine are outlined in the 


appendices. 
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i DISCUSSION OF RESULTS 


4.1 Single Cross-Tie Frames 


The results of the twenty-five frames evaluated 
in Appendix A varied significantly. The lightest frame 
had a cross-tie height of 12.0 meters (as did all of 
Lund's frames) and a side tank breadth of 12.0 meters. 
Lund's lightest frame had a side tank width of 16.0 
meters. The weights of the lightest limit design frame 
were 41.47 metric tons and 48.50 metric tons for factors 
of safety of 1.5 and 2.0 respectively. All of Lund's 
frames are based upon a factor of safety of 1.5 and his 
lightest structure weighed 42.1 metric tons. 

The heaviest limit design frame structure had 
a cross-tie height of 10.0 meters and a side tank width 
of 20.0 meters. The weights of this frame were 56.08 
metric tons and 63.45 metric tons for factors of safety 
of 1.5 and 2.0 respectively. 

The total frame weights of the elastic design 
structures and the limit design structures are within 
10.0 metric tons for all frames having factors of safety 
Oe li D4 

The distribution of frame weights resulting from 
the different design methods is significant. 

Member proportions were kept constant through 


both elastic and plastic limit design optimization, 


z= 


excluding plate flange areas. This procedure allows 
direct comparison of the relative weights of members 
with the same proportions by comparison of top flange 
areas. Comparisons of beam proportions are made using a 
factor Of safety of 1.5 for the limit design results 
unless stated to be otherwise. 

The limit design results show the size of 
Member #5, Figure A.l, to be larger than the Lund results 
for almost all frames considered. The limit design 
results with factors of safety of 1.5 gave top flange 
areas of 124 square centimeters in the lightest beam 
case and 283 square centimeters in the heaviest beam 
case. lLund's heaviest lower side tank beam had a top 
flange area of 127 square centimeters and his lightest 
had an area of 65 square centimeters. As might be 
expected the heaviest beams were those with the longest 
spans. The limit design lower side tank beams have 
roughly twice the cross sectional area of the elastic 
designs with the same side tank breadth in all cases. 
For example, Lund's frames with 16 meter side tank 
widths gave values of 92 and 127 square centimeters for 
top flange area. For similar frames the limit design 
mesults gave values of 198, 200, 203, 205, and 207 for 
top flange areas of the lower side tank beams. The 
dimensions of the Lund elastic frames may be found in the 


second rows of Tables 3 and 4 of Reference [5]. 
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The model used in the limit design study 
considered all loads acting on the structure. Lund 
does not discuss his model is detail. It is not known 
if he accounted for the influence of loads on the rigid 
extensions in his beam analysis. If loads in the rigid 
regions were neglected in the elastic designs, this 
could account for his smaller required size for the 
lower side tank beam. 

With one notable exception other beam elements 
were generally slightly lighter for the limit design 
case. The values calculated for the dimensions of 
Element #7 were much less in the limit design results. 
Representative values of from 69 to 101 square centimeters 
were Calculated for the top flange areas of the limit 
design frames versus a range of 134 to 253 for the 
elastic designs. In general, the elastic designs were 
roughly twice the weight of the plastic limit designs 
for this beam element (i.e. Element #7). This difference 
may have again been caused by unknown differences in 
models. 

All comparisons, so far, have been based on 
identical values of factors of safety of 1.5 for elastic 
and plastic limit designs. Figure A.5 shows that when 
a factor of safety of 2.0 is applied to the plastic limit 
designs, the total frame weights become higher than all 


but one of the frames designed by elastic methods. The 
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effect of increases in factors of safety on individual 
beam elements is illustrated in Figure A.4 for one beam 
element for the various plastic limit design frame 
geometries. Values of net frame area An are plotted 
against factors of safety for Member #1 of all frames 
evaluated. The dots represent values for the elastic 
designs of Lund. The increase in net area is seen to 

be nearly linear. On this figure only the line representing 
frames having cross-tie heights of 12.0 meters should be 
compared with the elastic design values since all elastic 
designs had 12.0 meter cross-tie heights. 

Plots for other beam elements would be similar 
to Figure A.4, but there could be up to 25 different 
lines if each of the frames considered had different 
values of Ap for the element considered. 

As would be expected from looking at Figure 
A.5, the weights of the individual limit design beam 
elements are in general heavier than Lund's results 
when a factor of safety of 2.0 is applied to the limit 
design results. 

The design effort required in evaluating the 
twenty-five frames discussed here was minimal. Total 
computation cost was $91.65. This results in a computation 
cost per frame of $3.66. The limit design problem was 
greatly simplified by using the assumption that the 


influence of brackets and beam overlap could be neglected. 
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Lund did not use this assumption. His elastic designs 
all assumed rigid joint regions. The elastic design 
computation time is much larger than for the plastic 
designs. Lund [5] cites CPU-time of 50 to 60 seconds 

on a UNIVAC 1108 for design of one frame structure. The 
limit design solutions of Appendix A required 
approximately 8 seconds in CPU-time for one frame 
structure on the IBM 370 including M.P.S.X. linear 
programming solutions, 

4.2 Zero Cross-Tie Frame with Rigid Corner Assumption 

The results obtained in Appendix B were for a 
Single frame based on the geometry of Lund's frame 
numbered OC22, 49.2t. The values of total frame weight 
for a number of factors of safety are plotted in Figure 
A.6. The weight of the elastic parent frame is plotted 
aS a point at a factor of safety of 1.5. From the figure 
it may be seen that the weight of the limit design frame 
calculated with a factor of safety of 1.5 is twenty per 
cent lighter than the elastic frame. Ata Factor of safety 
of 2.0, the limit design structure weighs approximately 
one half ton less than the elastic design structure. 

In the top line of Table 4, Reference Ee Lund 
lists the dimensions of his elastic design frame. Table 
A.2.4 lists values for the limit design frame for factors 
of safety of 1.0 to 2.5 in 0.5 increments. If the limit 


design dimensions with a factor of safety of 1.5 are 
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compared to the Lund results, it will be seen that the 
lower side tank girder, Member #2, is again substantially 
larger in the limit design results than in the elastic 
design results. Using top flange areas as a measure of 
relative weight, the limit design beam has a value of 

P32 cm’, while the elastic design beam has an area of 

167 cm*. This gives a relative difference in beam weights 
of approximately 15 per cent. 

As mentioned previously, it is not known if Lund 
accounted for the influence of loads on the rigid 
extensions. The limit design calculations accounted for 
all loads on the structure including the rigid regions. 

If Lund did not do this in his design model, his results 
could be lighter than the plastic limit design values. 

With the exception of Member #5, in the side tank 
top, all of the remaining limit design beam elements were 
about half the weight of their elastic design counterparts. 
Member #5 in the side tank top was about one fourth the 
weight of its elastic design counterpart. | 

The limit design calculations done in Appendix 
B were done by hand as opposed to Lund's computer approach. 
The hand operations were performed to allow close inspection 
of each mechanism. The results were taken to the computer 
for solution of the linear programming problem. 

Computation time was approximately 6 seconds. Again, the 


CPU-time for Lund's solutions averaged 50-60 seconds per 
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frame structure. The time required to do the limit 
design hand calculations was approximately 30 hours on 
the first time through. This time could be reduced to 
about four hours for successive runs by using the models 
already developed. 

These solutions could also be adapted for 
computer solution as was done in Appendix A using a model 
Similar to the one used in the hand calculations and a 
heavy reliance upon computer conditional tests. The 
conditional tests would be required to check the possible 


mechanisms for kinematic admissibility. 
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5. SUMMARY AND CONCLUSIONS 


An analytical study has been undertaken to 
investigate the use of the upper bound theorem of limit 
design in the design of transverse web frames in ships. 
Frames have been developed using limit design with 
results compared to frames of elastic design. 

In Appendix A the influence of beam overlap and 
brackets was neglected. This assumption was not made in 
the elastic design case. The results based on this 
assumption were frame structures of roughly equal 
weight to the elastic frame designs with some frames 
less than the least weight elastic designs if the same 
factor of safety is applied in both methods. For a factor 
of safety 2.0 the plastic limit design structures are 
approximately 10 per cent heavier. 

It is significant that while the overall 
structural weights are similar in this comparison, the 
distribution of weight is somewhat different. Frame 
structures designed by elastic methods have their 
largest members in one location while structures designed 
by the plastic method have their largest members in 
another location. 

In Appendix B the effect of beam overlap and 
bracket influence was idealized by assuming portions of 


these regions to be infinitely rigid. 
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Structural weight calculated using this assumption 
wasS approximately equal to the elastic design if a factor 
of safety of 2.0 was applied to the limit design solution, 
If a factor of safety of 1.5 is applied to the limit design 
solution, its result gives a twenty per cent reduction of 
we1ght Over the elastic design. 

A specific factor of safety was not specified 
for the plastic designs but rather a range of values was 
examined. The important advantage of limit design is that 
it allows the designer to define exactly the collapse 
moments for a given set of loading conditions. Factors 
of safety may be applied to these known quantities as 
required. 

In the case of the elastic design methods, 
factors of safety are defined somewhat arbitrarily 
without specific knowledge of when collapse will occur. 
The knowledge of specific collapse moments allows 
structures designed using plastic design methods to be 
much more efficient. 

The geometry of individual beam cross sections 
was not optimized or allowed to vary from the Lund 
proportions in this study. This was to permit direct 
comparison with the Lund designs. The limit design 
solutions could possibly have been made much lighter by 
use of more efficient beam cross section geometry. 


In the case where the assumption of rigid 
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corners was made, no attempt was made to calculate 
structural features within the rigid assumption zone. 

If this assumption is to be used, corner scantlings will 
have to be developed by separate analysis. 

Neglecting other possible advantages, the 
design costs uSing plastic limit design are very cheap. 
A design solution for a rather sophisticated tanker web 
frame can be had for very few dollars. This allows a 
number of different geometries to be examined as in 
Appendix A for the price of a single design found by 


elastic methods. 
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EIST OF FIGURES 


One bay portal frame. 


Fundamental collapse mechanisms of portal 
frame in Figure l. 


Tanker transverse web frame model with 
rigid extensions due to beam overlap and 
bracket influence. 


Corner model for calculation of rigid 
extensions due to beam overlap. 


Corner model with bracket for calculation 
of rigid extensions. 


Kinematically admissible collapse mechanism. 


Kinematically inadmissible collapse 
mechanism with w approximately equal to u. 


Graphical solution to linear programming 
problem of portal frame example. 


Tanker web frame model with one cross-tie. 


Loading conditions considered in tanker 
web frame design. 


Tanker web frame model without cross-ties 
and with rigid extensions. 


Frame cross sectional area vs. factor of 
safety for beam Element #1 of tanker web 
frame with a single cross-tie. 


Various frame weights for limit design and 
Lund elastic frames having a single cross- 
tie located twelve feet above the base line. 


Weights of zero cross-tie tanker web frame 


with rigid joints for various factors of 
safety. 
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Figure l. 


i ton 


One bay portal frame. 
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Figure 2. Fundamental collapse mechanisms 
of portal frame in Figure l. 
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Figure 3. Tanker transverse web frame model 
with rigid extensions due to beam 
overlap and bracket influence. 
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Figure 4a. Corner model for calculation 
of rigid extensions due to beam 
overlap. 
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Figure 4b. Corner model with bracket for 
calculation of rigid extensions. 
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Figure 5. Kinematically admissible collapse 
mechanism. 





Figure 6. Kinematically inadmissible 
collapse mechanism with w 
approximately equal to u. 
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Figure 7. Graphical solution to linear 
programming problem of portal 
frame example. 
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Figure A.l. Tanker web frame model with 
one cross-tie. 
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Load #1 Load #2 





Load #3 Load #4 


(Loads in kp/cm) 


Figure A.2. Loading conditions considered in tanker 
web frame design. 
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Figure A.3. Tanker web frame model without 
cross-ties and with rigid 
extensions. 
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Figure A.4 Frame cross sectional area vs. factor of 
safety for beam Element #1 of tanker web 
frame with a single cross-tie. 
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Figure A.5 Various frame weights for limit design and 
Lund elastic frames having a single cross- 
tie located twelve feet above the base line. 
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Figure A.6 Weights of zero cross-tie tanker web frame 
with rigid joints for various factors of 
safety. 
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APPENDIX A - TANKER FRAME COMPARISON 


Modern supertankers are designed with large 
variations from ship to ship in the internal primary 
stiffening system. The number of cross-tieS varies from 
zero to three. Structures may be designed with or without 
longitudinal side girders. Bottom center girders may 
be low or high. When compared on the basis of weight 
to strength ratios, some designs should be superior to 
others [5]. 

The present study comprises investigations 
of variations of a single transverse web frame 
configuration. 

A FORTRAN computer program was developed to 
evaluate the limit moment values of a tanker webframe 
having a single cross-tie located in the wing tank. 

The vessel considered is identical to the one considered 
by Lund in Reference [5]. The vessel is a typical 


supertanker with the following principal dimensions: 


Length Lop = 320.0 m 
Breadth B= 50.0m 
Depth D= 25.0m 
Draft qd= 19.5 m 
Camber (all in side tank) = 1.0m 
Spacing of transverse frames = 5.0: 1m 


The basic geometry of the frame studied is 
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shown in Figure A.l. The program was developed 
specifically to evaluate variations of this frame, but 
it could be made more general without extenSive 
modification. As written, the required input values 
are the cross-tie height above the base line S2 and 
the wing tank breadth S5, both values being in 
centimeters. All other values are either written 
into or calculated by the program as needed, based 

on the parent vessel geometry and the four loading 
conditions prescribed by Lund [5] and to be described 
later in this appendix. System lines are assumed to 
be formed by the inside surface of the shell plating 
and the center tank face plating of the longitudinal 
bulkhead. It was found when beams were sized that 
plastic neutral axes for the various beam elements 
remained at or very near these system lines. 

The possible influence of beam overlap and 
brackets was neglected with plastic hinges allowed to 
form anywhere along the lengths of individual beam 
elements. Combined plastic collapse mechanisms which are 
very complex and are unlikely to be actual collapse 
mechanisms have been eliminated from computation 
a priori to reduce computing costs. 

The range of cross-tie heights examined is 
from 10 meters to 14 meters in one meter intervals. 


The range of side tank breadths examined is from 12 
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meters to 20 meters in 2 meter intervals. In all a total 
of twenty-five separate variations of the same single 
cross-tie web frame were examined. 

The loading conditions applied were identical 
to those used by Lund in his elastic designs. The 
various conditions are shown in Figure A.2. The loading 
models as applied represent the following situations 
54. 

1. Empty center tank, cargo in side tanks, full draft 
with dynamic additions. 

2. Empty side tanks, cargo in center tank, full draft 
with dynamic additions. 

3. Empty side tanks, ballast in center tank, draft 

Of 70.25. Ds 

4, Ballast in center and side tanks, draft of 0.35 D. 

As stated previously, the program input 
variables are cross-tie height and side tank breadth, 
S2 and S5 respectively. The input format is free 
with the real values of S2 and S5 in centimeters listed 
in the same order shown here. The program sequence 
exits when a zero value for S2 is read so the last 
data card should have real values of zero and any other 
real constant listed consecutively. 

The output values of the FORTRAN program 
include the coefficients of the internal energy 


dissipation portion of the limit design upper bound 
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theorem. These coefficients are presented as the 

values as 3 Gb Qthe scOeLrIclent Mater Ain the limit 

design problem constraint relation. The associated 

values of the external work rates in the principle of 
virtual work relations are also given as output. These 
work values represent the right hand side of the constraint 
equations b;- 

The output information is all punched on 
cards to be used as input data in the linear programming 
problem to be solved by the M.P.S.X. subroutine. 

The format of the punched cards is quite important, 
but is not discussed here. A discussion of M.P.S.X. 
input format is given in the numerical routine appendix. 

Finally, a punched card is inserted at the 
front of each output data deck listing the values of 
S2 and S5. This is to enable the user to identify 
the various output data sets when more than one frame 
is examined. 

The program was checked for logical and 
numerical errors by comparing output values and hand 
calculation results on the web frame having S2 = 12.0 
meters and S5 = 12.0 meters. 

The complete listing of the program with 
sample data values is presented on the pages that 


follow. 
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"T= (2‘7ZT)Wv 
*p/Odx (Z¥*SS=(TT)M 
SV = eS iy 
*STT=0d 
(ZJx¥%6S) xdd=(O0T)M 
"p=(6°0OT) WY 
*S80T=dd 

AONENTANI T GWOT O 
8S=6S 
GS-°00Sz=8S 
G°xx((Z¥%SGS)+°T)=€S 


SD=Vo 
LS-°009Z2=9S 
CS=LS 
cS-°00SZ=TS 


NOLLVOIVARZ HLONAI LNAWAIG O 

C OLOD(°0°Os° ZS) AL 

cs‘zs‘awaua T 

(66) M‘ (6°66) VY NOISNAWIG 
qors 
OO0OO0Z=HZISM1d=d0d‘s ad NISAS*O// 
AIALWM OSX // 
(Q4EMMNWL‘S ‘4 8067Td‘ 9TSETN) =YRSN GILINs/ 
MOSZ=NO1Dgua ‘ V=SSWIO’ ,NHONIM, // 


T°w eaznbtg 
UT TSPpOW sy UO paseg sWeAA qom AsyueL, Fo ubtsseq 4Atwty AoFJ werzborg NWULuow 


-6§4— 





°9/( (Wdx° @+dd) «Zxx% (HS-TS) +(HS-TS) /HS«* (Xd+Wds°Z) ¥xZx*HS)=(LI)M 
Xd+(HS-TS) /(Xd-dd) »HS=Wd 
((HS-TS) /HS) ¥°Z+°Z=(T‘ LI) ¥ 

(Wx °Z) /(G° xx (Oxdx * 0-7 x0) -d-) =HS 
*ZT/(Xdx° Z+dd) x TS=0 
°9/(ddx°7+Xd) -°€/ (TS) =a 
°7/a-(TSx%°b) /(Xd-dd) =4 
(ST)M=(9T)M 

"T=(€‘9T) 

(ZS/HS) «°Z+1TT=(2‘9T) © 

°9/( (HS—ZS) /HSx (Wdx °Z7+Xd) xZxx (HS—-ZS) + (Hd+Wdx°Z) «Z%*HS) =(ST)M 
Hd+7dS/ (Hd-Xd) xHS=Wd 
((HS-ZS) /HS) «*°Z+°7Z=(7Z‘ST)W 
(Wx°Z)/(G° xx (Oxdx *° ¥-Z4%d) -G-) =HS 
*ZL/ (Hdx °Z+Xd) x«ZS=O 
°9/(Xdx°Z+Hd) -°€/(ZSxUu) =a 

(°7/a) -(2S%°t) / (Hd-Xd) =4 
(TS+ZS) / (Hd-Xd) =u 
Hd+(ZS+TS) /( (Hd-dd) *ZS) =Xd 
‘OVI=dd 

*STI=Hd 

(€T)M=(PT)M 

*7=(6‘ PI) WY 

iL. a) a 

"€=(G‘PT)Y 

(TT) M+°Z/ (OT) M=(ET)M 

“T=(6ET) WY 

o1—(f ew 

-2=(G S10) w 

1l—(2 a) a 

(IT)M= (ZI)M 

= Gch) 
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(W¥°Z)/ (Go x (Ox ge ¥-Z¥ 4%) -G-=HS 
*ZT/ (Mds ° Z+dd) ¥LS=0 

°9/ (dds °Z+Md) -°€/(LS%u) =a 
°7/a-(LS%°P) / (Md-dd) =4 
°9/((HS-9S) / (Wd *Z7+Md) ¥HS¥Z*% (HS-9S) +(Ad+Wdy°Z) «Zx¥HS) =(TZ)M 
(HS-9S) /HS¥°Z+°Z=(9‘/TZ)W 
Ad+9S/(Ad=-Md) x»HS=Wd 
(W¥°Z)/(S°)) (Ox de ° P-Ze¥%G) -G-) =HS 
*ZT/ (Adx ° Z+Md) ¥9S=0 

*9/ (Mdx *Z+Ad) -°€/(9S%u) =4 
°7/a-(9S%°) /(Ad-Md) =4 

(LS+9S) /(Ad-dd) =u 
Ad+(LS+9S) / (Ad-dd) ¥9S=Md 

007 T=dd 

*S6=Ad 

(6T)M=(0Z)M 

(HS-€S) /HS+°Z=(€/0Z) WV 

(HS-€S) /HS=(T‘0Z) W 

°9/( (HS-€S) /HSx (Ldx° Z+Wd) ¥Z¥% (HS-€S) +(Ld+Wdy °Z) xZ¥xHS)=(61)M 
(HS-€S) /HS¥°7+°Z=(€’6T) WY 
Ld+(Ld-Od) x (€S/HS) =Wd 

(W¥°Z)/ (So xx (Ox de ° ¥-Z¥%d) -A-) =HS 
*ZT/ (Ld °Z+0d) x€S=0 
°9/(Odx °7+Ld) -°€/ES x= 
°7/a-(€Sx%° Pb) /(Ld-Od) =4 
€S/(Ld-Od) =u 

*S¢6=Ld 

“OVT=0d 

(LT) M=(8T)M 

*T=(€’°8T)W 

(HS-TS) /HSx¥°Z+°T=(T’ 8T) ¥ 
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(ZS+1TS) /Hd=u 
(ZS+TS) /Hdx TS=Xd 
*S$80T=Hd 
(LZ)M=(8Z) M 
"T=(6'82) W 
EGE) rd 
"€=(G‘82) 0 
°Z/ (97) M+ (GZ) M=(LZ) M 
wa B= Sa brah ag 
“T=(L‘LZ) W 
"T=(G‘LZ)W 
“Ute 3c) © 
(SZ) M=(97) M 
"T=(7Z‘972) © 
"€=(G¢‘97Z) WV 
°P/Tx*xGSx0d=(GZ)M 
"p=(G‘SzZ)W 
*G80T=0d 
"0/7" (6S%°7) xVd=(972)M 
"p=(6'97) WY 
*S{TT=Wd 

GAONANIANI Z AWOT O 
°9/( (Xd °Z-Ld) xZxxTS-LS/T 
9Sx% (Hd+Xdx°Z) xZ¥%x7S— (Ld+Mdx ° 7) xZ¥%9S+LS/9S% (Mdx °7-dd) ¥Zxx%1S) =(E€Z)M 
LS/9S=(L‘E7) WT 
CS/TS=(7‘4 ET) 
LS/9Sx°Z=(9'E7) W 
LS/9Sx%°7=(PE7) 
LS/9S+°7=(T‘4EZ) WY 
°9((HS-LS) /HS»% (Ndx °Z+dd) «Zx% (HS-LS) + (Md+Wdx°Z) x °Zx¥xHS) =(ZZ)M 
(HS-LS) /HSx°Z+°7=(L‘477) W 
Md+LS/ (Md—-dd) »HS=Wd 
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(LS+9S) / (Ad-Wd) =4 

“007 T=Wd 

°c6=Ad 

*p/Zxx%(8S%°7Z) xAd=(SE)M 

“p=(8‘SGE) WY 

°G6=Ad 

(€£)M=(P7E)M 

“T=(€°pE)Y 

HS/(HS-TS) «°Z+° T=(T’ PE) WY 

°O/ (Hdx* ZxZxx (HS-TS) +HS/ (Wdx °7+Xd) » (HS-TS) «Z¥%HS) =(€£)M 
TS/ (HS-TS) Xd=Wd 

HS/(HS-TS) «°Z+°Z=(T‘ EE) W 

(G*xx(°€/°T) -°T) x TS=HS 

SCRE 

/ (Xd+Hdx"Z) xZS¥SSt(7Sx%°) /GSaTax (°7/GS) xOd+* Z/Zx%(°7Z/GS) xDd=(ZE)M 
7S/GS+°7Z=(S‘ZE) WY 

(ZS%°Z) /SS=(2‘ZE) W 

°9(SS/HSZ 

x (Wdx° 7+Hd) «Z%%HS-ZS) +(Xd+Wdx°Z) «Ze¥%HS) +(°7Z*%SS) /HS¥Z¥xGSxOd=(TE)M 
GS/HS=(S‘TE)W¥ 

GS/HS+°Z=(Z‘TE) 

(6Z)M=(0€)M 

“T=(G‘O0E) WV 

HS/ (HS-ZS) *°Z+° T=(Z‘0€) © 

°9/( (Wd ° Z+Hd) xZ*xHS-ZS) +HS/ (HS—-ZS) » (Xd+Wdx°Z) *«Z¥xHS) =(6Z)M 
Xd+2S/ (Xd-Hd) xHS=Wd 

HS/(HS-ZS) «°Z+°Z=(2‘6Z) W 

(Wx °Z)/(G* xx (Dx de ° $-Z¥%0) -A-) =HS 

“ZL/(Xdx%°Z+Hd) ¥ZS=09 

°9/ (Hdx °Z+Xd) -°€/72Sxu=a 

°7/a-(7S%°f) /(Xd-Hd) =4 
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“T=(9¢‘Zp) Ww 
°9/Hd*Z**0=(TP) M 
(HS-ZS) /HS¥°Z+°Z=(Z‘ TP) W 
(S*xx(ZS/d) -°T) =HS 
°OZE=Hd 
°Sz9=d 
(6£)M=(07)M 
"€=(9'07) ¥ 
“T=(2°07) W 
*p/OduZ4%%GS=(6£)M 
“p=(G‘6E) 0 
“O7E=0d 
*p/WduZe%(6S%°Z)=(8E)M 
*cTOT=Wd 
“p=(6°8E) WY 
AONANTANI € GWOT DO 
°9/((HS-LS) / (Wd«°Z+Wd) «HS*Ze%(HS-LS) +(Md+Wdx°Z) *Z*%*%HS) =(LE)M 
(HS-ZLS) /HS#°Z+°Z=(L‘ LE) © 
Md+LS/ (Md-Wd) HS=Wd 
(Wx¥°Z)/ (Go ex (Dede °P-ZJ¥%d) -d-) =HS 
*ZT/ (Md °Z+Wd) xLS=O 
°QO/ (Wd °Z+Md) -°€/LSxu=d 
°7/a-(LS%°%) 7 (Md-Wd) =4 
°9/( (HS-9S) /(Wdx °7+Md) +HS*Z x (HS—-9S) + (Ad+Wdx°Z) *Z¥*HS) =(9€)M 
(HS-9S) /HS*°Z+°Z=(9‘/9€) 0 
Ad+9S/(Ad-Md) xHS=Wd 
(Wx °Z)/ (Go x (Oxde° P-C4¥%) -d-) =HS 
“ZT/ (Ads °Z+Md) ¥9S=9 
°9/ (Md °Z+Ad) -°€/9S#u=4 
°¢/a-(9S%°P) / (Ad-Md) =J 
Adt+tidx9S=Md 
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(HS-TS) /HS%°Z+°Z=(2‘ 8h) 0 
Wd+(°S/8-72S+1TS/ (Od-Dd) =Xd 
°0S=0d 
°GL8=9d 
TS*€S °O=HS 
(97) M=(LU)M 
"€=(G‘LP)W 
“T=(Z‘Ld) W 
°P/Tx%SGSxdd=(97)M 
“p=(S'9P) 
°¢L8=dd 
"p/Tux(6S%°Z) xWd=(SP)M 
*p=(6'SP) 
°¢L8=Wd 

AONANTANI ~ GWOT O 
°9/((HS-LS) /(Wdx* Z+Wd) +HS%Z%% (HS-LS) +(Md+Wdx°Z) xZ¥*HS)=(P7P)M 
(HS-LS) /HS%°@+°Z=(L‘ PP) © 
(HS+9S) xd¥=Wd 
(Wx °Z) /(G° ux (Oxde° 7-Zx%d) -G—-) =HS 
*ZT/ (Md °Z+Wd) x LS=O 
°Q9/ (Wdx °Z+Md) -°€/LSy¥u=a 
°Z/a-(LS%° 7) /(Md-Wd) =4 
(L£S+9S) /Wd=u 
°9/((HS-9S) /(Wdx°Z+Md) x€%xHS+Wd¥Z xx (HS-9S) )=(€7)M 
(HS-9S) /HS¥°Z+°Z=(9‘EP) © 
(LS+9S) /(HS-9S) xWd=Wd 
(G*xx(°€/°T)-°T) «9S=HS 
(2S+9S) /9S*xVd=Md 
*SEETHVd 
*0=Ad 
(Tp) M=(ZP)M 
(HS-ZS) /HS¥°Z+°T=(7‘4Z7) 0 


-70-= 





(Ri) ye “OTP noNna 

G OL OD (°O°OR* UAWNW) ALI 

(c’I) W=aaWnw 

Ce Shh bea Od 

6°T=£ 9 od 

GS‘zs‘ZT HONNd 
*7/dd¥Z%xGS=(€S)M 

"Z=(S‘ES) WV 

"f=(P°ES) WY 

*Z=(€°ES) WV 

*S80T=dd 

°9/ (Wd °Zx (HS-€S) + (Wdx °Z+Od) ¥+Z¥¥HS) =(ZS)M 
€S/(HS-E€S) xOd=wWd 

(HS-€S) /HS¥°7+°Z=(€'4ZS) WV 
(S*xx(°E/°T) -°T) x€S=HS 

°O0=Ld 

°0S=0d 

(0S)M=(TS)M 

"T—(Ce. 1S) 2 

(HS-TS) /HS%°2+°T=(1T‘TS) 
°9/((HS-TS) / (Ndx °Z+Xd) ¥Z¥% (HS-TS) xHS+ (Od+Wdy°Z) ¥Zx*%HS)=(0S)M 
(HS-TS) /HS¥°2+°Z2=(T‘0S) V7 

Od+HS #¥=Wd 

(Wx °Z)/ (Go 4% (Dx dx °7-Zx%¥d) -d-) =HS 
*ZT/(Odx °Z+Xd) x TS=9 

°9/ (Xdx°Z+0Od) -°€/TSxu=4 
°Z/a-(TSx%°) /(O0d-Xd) =a 

°SZ9OT/ (Od-Dd) =u 

(87) M=(69)M 

(HS-z2S) /HS=(S‘6%) W 
(HS-ZS) /HS+° Z=(7‘6h) © 

°2/Ddx (HS-TS) +°Z/Ddx¥Zx¥xHS=(87)M 
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The output data from the program is not listed 
in this document as computed. To list this information 
would require a separate volume. 

To complete the M.P.S.X. linear programming 
data set, objective function cards listing the various 
beam lengths are added to the FORTRAN program output 
cards. Finally the desired M.P.S.X. control language 
1s inserted and the problem run to get optimum upper 
bound values for the plastic limit moment for each 
beam in each frame. It should again be noted that the 
very rigid M.P.S.X. input format requirements outlined 
in the numerical routine appendix must be followed 
precisely. Failure to do so may result in the linear 
programming routine being abandoned or worse, in 
erroneous results. 

In order to compare the web frame sizes required 
by the plastic limit design problem to Lund's elastic 
design results, the values of plastic limit moment are 
used to size beam elements based on the beam proportions 
used by Lund in his analysis. A 0.2 cm corrosion 
allowance is then applied to plating thicknesses. Lund 
based his beam proportions on a set of starting values 
which he scaled up to meet stress control requirements. 
The values of plate flange area a are likely to be 
drived by design factors not considered here and are 


assumed to be held constant. Lund does not list thickness 


he 





values for plate flanges. It is therefore assumed that 
the full area of A, may be considered to be part of 
the beam to be sized. The starting values of web 
depth, web thickness, top flange area and plate flange 
area as used by Lund and in this study for scaling beams 
is listed in Table A.1l.l. 

Lund assumed his designs to be of mild steel. 
Stress values were restricted to 1500 kp/em? versus the 
yield stress value of 2249.8 kp/cm*. This gives an 
effective factor of safety of 1.50. There is some 
question as to the "proper" factor of safety to be 
applied to the values of plastic limit moment to make 
a reasonable comparison between the elastic and plastic 
designs. 

Values of beam geometry for the plastic limit 
designs were calculated using factors of safety of 1.5, 
2.0, and 2.5. The dimensions calculated included those 
listed in the heading of Table A.1.1, as well as the 
total cross sectional area of the beam excluding plate 
flange area. This total cross sectional area will 
henceforth be called Ane 

In calculating the values named above the 
following relationships were satisfied. 

M, = 0. 2 


where 
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In these relationships M, is the plastic limit 
moment calculated in the limit design problem, ws is 
the material yield stress, A is the total required 
area of the cross section and a is the distance between 
centroids of the two equal areas lying above and below 
the plastic neutral axis. 

The values Ans web depth h, web thickness t 


and top flange area A are listed in the tables at the 


? 
end of this appendix, starting with Table A.1.2. The 
numbers in the frame column represent values of S2 

and S5 respectively. For example "12 - 14" denotes that 
this is the frame where S2 = 12.0 meters and S5 = 14.0 
meters. Dimensions with factors of safety of 1.5 and 
2.0 are listed for all beam elements exclusive of 
Element #4. It was found that the limiting mechanisms 
of the limit design problem never included this element. 
Those mechanisms in which failure of the member might 
have occurred had high values of internal energy 
dissipation with low external work rates when compared 
with the other possible collapse mechanisms. 

It is known that Member #4 is in compression 
and subject to buckling. Plastic buckling is not 
considered in this study. In order to provide some 
measure of comparison, the dimensions found by Lund 
in his study are used unaltered in this analysis when 


final frame weights are computed. In these computations 
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Element #4 is assumed to have a cross sectional area of 
432 cm? with a web height of 153 cm for all cases, 

Figure A.4 illustrates variations of net 
frame area An for various factors of safety as calculated 
for Member #1. 

Total web frame weights were calculated for 
one half of the complete frames. To be consistent with 
Lund, comparisons are based on the assumption that the 
relative amount of local stiffening and relative bracket 
sizes are equal for all investigated frames and these 
components are neglected. Overlaps are compensated 
for by subtracting one half of the average weight in 
the joint region for each beam. Total web frame weights 
were therefore calculated using the following relation: 

9 
W = p) {1, A, - % (ec; + c,) Agh 
i=1 
where p 1s material density, 1; 1s the member length 
and C5 and Cc, are the overlaps at the respective beam 
ends. This relationship was also used by Lund [5]. 

The weight in metric tons of each frame 
considered is listed in Table A.1.18. Weights were 
calculated based on factors of safety of 1.5 and 2.0. 

Lund assumed a constant cross-tie height of 
12.0 meters. The limit design frame group with the 
12.0 meter cross-tie height had the lowest weights, 


aS a group. Figure A.5 shows values of frame weight 
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versus side tank breadth for single cross-tie frames 
with cross-tie heights of 12.0 meters. Solid lines 
represent limit design results. Dashed lines represent 
all values for single cross-tie frames given by Lund. 
There are two curves joining the elastic frame weights. 
The one to the upper right represents weight values 
where the bulkhead girders are outboard of the 
longitudinal bulkhead. The remaining elastic design 
curve gives weight values for frames with bulkhead 
girders inboard of the longitudinal bulkhead. The two 
limit design curves both assume the bulkhead girders to 
to inboard of the longitudinal bulkhead. The two 
curves give weights for factors of safety of 1.5 and 


2.0 for the various frames. 


ec 





TABLE A.1.1 


Initial Member Dimensions on which 
Beam Proportions are Based 


Web Top Flange Plate Flange 

Member Height Thickness Area Area 
No. (cm) (cm) (cm2) (cm) 
1 300 leo 120 660 
2 300 1S: 120 660 
3 300 14 80 750 
4 160 2.0 150 150 
5 400 1.8 120 750 
6 300 1S, 120 420 
a. 300 lee, 120 510 
8 240 15. 80 7350 


9 400 2.0 120 Fesyo 
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Frame 


(S2 = S5) 


10-12 
10-14 
10-16 
10-18 
10-20 
11-12 
11-14 
11-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
13-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


(Factor Safety 


A 


T 
(cm) 


Og 


is Ko a 


Table A.1.2 - Dimensions of Member 
125) 


#1 





Frame 


(S2 - S5) 


10-12 
10-14 
10-16 
10-18 
10-20 
L1l-12 
11-14 
11-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
13-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


Table A.1.3 - Dimensions of Member #2 


(Factor Safety = 1.5) 


A 
(cm2) 


e130 
642.58 
669.92 
699.33 
726 —a) © 
658.35 
687.79 
(lse22 
748.63 
782.15 
701.43 
733.95 
766.44 
UR els 
833.42 
745.49 
780.06 
BS S50 
848.05 
884.61 
771.68 
811.46 
B59.55 
897.14 
932.61 


h 


(cm) 


sve oR Sey 
302.88 
309.58 
316.64 
3235.55 
306.77 
513.09 
321.10 
328.15 
339,16 
317.14 
Pe Sa) 
BSiace2 
339.74 
347.09 
327.43 
B35, 9 
342.74 
350.26 
358.06 
333.40 
342.28 
392573 
360.70 
368.07 


=op— 


(cm) 


1.68 
eye 
uf i) 
NTs: 
Eo 
gee | 
LUT AF 
io. 
1.84 
Neck: 
L419 
eo 
1.86 
UBS 6) 
Be 
1.84 
Too 
iol 


ie 


I) 
1.87 
SOUL 
1.96 
2.00 
2.04 


A 
(cm2) 


116.42 
Pees 3¢ 
4 eho 
133.68 
139255 
425.47 
DBs es 
137.47 
143.58 
5 U2 
134.41 
140.63 
147.16 
133.39 
160.63 
142.95 
149.09 
156.63 
103.56 
LO. 95 
148.21 
LD O6 aL 
165.69 
173.47 
180.63 





Table A.1.4 -—- Dimensions of Member #3 
(Factor Safety = 1.5) 


Frame Am h sof Ap 
(S2 - S5) (cm?) (cm) (cm) (cm) 





10-12 426.89 259.46 1.41 59.84 
10-14 " " 

10-16 7) " " " 
10-18 " " " " 
10-20 : " Ww " 
11-12 386.19 245.93 1.35 53.76 
11-14 " 

11-16 " " '" 
11-18 1" 1" 1 " 
11-20 ” " " oe 
12-12 346.40 231.99 1.28 47.84 
12-14 n " 

12=1G 1 1 " 1" 
15-18 iT Ww " " 
12-00 1 ut ut " 
13-12 309.68 218.40 22 42.40 
13-14 q 
13=16 1" i) " . " 
13-18 iT) " 1"! " 
1o=20 " " 1" " 
14-12 113.03 125.14 ~ 78 VES Eps 
14-14 35.76 1385.7 8 < 85 Ve ek2 
14-16 273.07 204.35 Zao S512 
14-18 " 7) 
14-20 1" 1" " " 


=e] = 





Frame 


(S2 =- S5) 


10-12 
10-14 
10-16 
10-18 
10-20 
11-12 
11-14 
11-16 
11-18 
11=20 
12-12 
12-14 
12-16 
12-18 
12-20 
13-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


A 


(cm2) 


956.60 
V21665.0 
1488.72 
L766.222 
2048.14 

980.67 
L237 532 
1508.43 
1783.80 
2064.64 
1003.68 
EZ LAOS 
1526.05 
1800.34 
2079.08 
1022.49 
1274.78 
1541.60 
1813.78 
2092.48 
LOGS. 77 
1338.20 
lhe sy oye al 
1827.21 
2104.85 


(cm) 


408.01 
462.50 
SS 02 
Dole te 
605.69 
413.35 
466.60 
Sl? lS 
5964.00 
608.20 
418.39 
470.46 
520.26 
566.69 
610.39 
422.46 
473.89 
923.00 
568.88 
612.41 
438.83 
485.99 
525,594 
7 £0 
614.27 


ER oy a 


Table A.1.5 - Dimensions of Member #5 


(Factor Safety iD) 


(cm) 





2.04 
2.20 
2 
Vere 
25.95 
2.06 
2.50 
2553 
2.74 
2294 
2.08 
2552 
2.54 
Lie PD 
PP o/s 
20 
Pe 
2.52 


Ze7O 


2.96 
Zenley 
Zao 
Za50 
Paid 
PAP EN 


A 
(cm2) 


124.86 
160.43 
IES eae eie: 
236.14 
ZIS<¢i 4 
2 ee 
63.29 
200557 
238.57 
271,43 
dN B74 
166.00 
203.900 
240.86 
2/9545 
133.86 
168.43 
205.14 
24207 | 
281.29 
144,43 
77 2 ee 
207.14 
244.57 
233300 





Frame 


(S2 - S5) 


10-12 
10-14 
10-16 
10-18 
10=—20 
ll-12 
11-14 
11-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
13-12 
13-14 
13-16 
13=-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


Table A.1.6 - Dimensions of Member 


(Factor Safety = 1.5) 


(om2) (em) 

384.13 230.67 
ih ia 
a a 
aa re 


234.36 PyG,.6 1 


eee 








Table A.1.7 -— Dimensions of Member #7 
(Factor Safety = 1.5) 


Frame Am h te Ap 
(S2 —- S5) (cm2) (cm) (cm) (cm2) 


LQ=-12 372.56 226.60 oS 68.63 


10-14 
10-16 "W " " 
10-18 " " " 
10-20 " " " 


11-12 414.01 240.07 1.40 76.84 


11-14 
tledke " " "W 
11-18 " " " 
1LP=90 1" " " 
2 EZ 454.45 Z52e22 1.46 84 
12-14 tt Ww 
19-16 "W " " 
12-18 "W " " 
12-20 " " "W 
ke hy 494.78 263.88 I Sy 92 
Ikea " " 
13216 "W " " 
1 3=18 "W " " 
13-20 " " " 
14-12 5935.00 275.01 1.58 100 
14-14 " " 
14-16 " " " 
14-18 " " " 
14-20 " " " 
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Table A.1.8 - Dimensions of Member #8 
(Factor Safety = 1.5) 





Frame Ap h ac Ay 
(SJ. SS) (cm2 ) (cm) (cm) (cm2) 

10-12 27267 188.33 1.18 49.26 
10-14 221.67 168.36 1.08 39.37 
10-16 173.45 147.23 GT efoealal 
10-18 127 207 124,34 285 AF 
10-20 85.85 99.26 a72 13.68 
ae 272267 188.33 1.18 49.26 
11-14 221.67 168.36 1.08 39.37 
11-16 173.45 147.23 .97 30-l 
11-18 27 287 124.34 £85 DUS AF 
11-20 85.85 99.26 a, 13.68 
12-12 272.67 188.33 1.18 49.26 
12-14 221.67 168. 36 1.08 39.37 
12-16 173.45 147.23 .97 Sone) 
12-18 Te 124.34 £85 2 Ag 
12-20 85.85 99.26 Te 13.68 
1312 Doo 188.33 iL IL 49.26 
13-14 DOI GF. 168. 36 1.08 39.37 
13-16 173.45 147.23 .97 301 
13-18 127-97 124.34 285 Blea y 
13-20 85.85 99.26 ae 13.68 
Ae? DO AGE 188. 33 1.18 49.26 
14-14 221.67 168.36 1.08 39.37 
14-16 173.45 14g 222 .97 g021 
14-18 127.87 124.34 285 21.47 
14-20 85.85 99.26 72 13.68 
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Frame 


(S2 - S5) 


10-12 
10-14 
10-16 
10-18 
10-20 
11-12 
11-14 
11-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
13-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


(Factor Safety 


A 


(cin2) 


1278.87 
Osi. 25> 
800.16 
581.14 
381.06 
EZ7 C607 
L031. 25 
800.16 
581.14 
381.06 
2 Cao) 
OSts. 25 
800.16 
581.14 
381.06 
VZ278.07 
Oost .25 
800.16 
581.14 
381.06 
Zao 
1031.25 
800.16 
581.14 
381.06 


h 


(cm) 


454525 
406.26 
350.02 
S00. 72 
240.29 
454.35 
406.26 
SS 3) 
B00. 2 
240.29 
454.35 
406.26 
Sanaa S74 
500.9 2 
240, 2.0 
454.35 
406.26 
3 Oo Se 0 2 
300272 
240.29 
454,35 
406.26 
SSP 874 
300.72 
240.29 
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Table A.1.9 - Dimensions of Member #9 
5) 


(cm) 


Zea) 
2623 
eG 
ea v 
1.40 
Aa: 
2.25 
1.98 
ike 0 
1.40 
2.47 
Para 
aoe 
Pea 9 
1.40 
2.47 
2,25 
1.98 
1.70 
1.40 
2.47 
Zeca 
1.98 
ies 0 
1.40 


A 
(om?) 


£547. 63 
25-0 0 
94:95 
67283 
43.30 
Leas 
2316 
94.96 
6/285 
43.30 
154.83 
23476 
94.96 
67.83 
oS U 
154-05 
1235 
94.96 
67.83 
43.30 
154.83 
ZS iene 
94.96 
67.83 
A300 





Frame 


(S2 = S5) 


LO= 12 
10-14 
10-16 
Shea. 
10-20 
11-12 
11-14 
11-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
3=12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


Table A.1.10 -—- Dimensions of Member #1 
(Factor Safety = 2.0) 


Ap h t 

(cm2) (cm) (cm) 
514.90 269.50 eas > 
465.07 255.36 1.48 
417.21 24-205 | ae ae 
372.38 226.88 le33 
370624 226.18 e353 
364.89 224.43 e352 
329750 212250 i. Z26 


i a 


86.99 





Frame 


(S2 - S5) 


10-12 
10-14 
10-16 
10-18 
10-20 
11-12 
11-14 
11-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
3-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


(Factor Safety 


ia 
(cm2) 


7307.05 
770.63 
804.13 
840.74 
876.26 
790753 
826.10 
862.68 
900.27 
38971 
842.83 
881.48 
921.14 
960.76 
LOO 32 
895.05 
936.78 
977.44 
1020.13 
1038.87 
927.40 
OF S32 
1033.66 
10738, 39 
1122.05 


(cm) 


S22.) 
33357 
340.66 
348.68 
356.30 
337.64 
345.50 
353.40 
361.36 
SiG). 57 
549. 13 
357.40 
365.70 
373 ou: 
Bie. © 
360.26 
BGG. 02 
OT 8 
385.66 
3094.53 
366.99 
376.76 
388.32 
396.96 
405.23 
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Table A.1.11 - Dimensions of Member #2 
2.0) 


(cm) 





1s 83 
7 
1.90 
1.94 
9s 
189 
do 
Pero 
22.01 
2.00 
US), 


208 
Zu 
Zeer 
2500 
2.04 


PAS ORS ae 


Zo 
Zi 
2203 
2.08 
2.14 
2.18 
223 


A 
fame) 


cad, o> 
148.00 
154.74 
aE Zielab 
UGDLZ6 
152.100 
159). 6 
166.53 
Oley a: aro | 
G2 jae 
L622 33 
170.32 
1 Be332 
66.32 
94.53 
7 30D 
181.47 
69.08 
eS Fe eb 
2022 11 
1795016 
89.26 
201 20> 
220 ae 
PAI > oe) 





Frame 


(S29 = S5) 


10-12 
10-14 
10-16 
10-18 
10-20 
11-12 
11-14 
11-16 
11-18 
Li=-20 
12-12 
12-14 
12-16 
12-18 
12-20 
L3=-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


Table A.1.12 - Dimensions of Member #3 
(Factor Safety = 2.0) 





Aq h t 
(cm?) (cm) (cm) 
Siva 235255 Le53 
462.13 270.67 1.46 
i i we 
A2USe 2 255-62 1.39 
tf Ww te 
it WwW 7 
370.07 240,37 ae 
1? uw wT 
uf ia uw 
134-63 Bey see es: 34 
161.59 L522 297) 91 
326.98 224.90 a 25 
1 Ww Tt 
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58.08 


Die ob 





Frame 


(S2 - S5) 


BO=12 
10-14 
10-16 
10-18 
10-20 
11-12 
11-14 
Fl=-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
P3=-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


(Factor Safety 


A 
(cmn2) 


1149.83 
1462.79 
VIZ 07 
2126.50 
2467.32 
1133.80 
1488.72 
1815.84 
2148.15 
2486.87 
1206.09 
LoL2Z 757 
1837.54 
ZG 77S 
2505.58 
1230.03 
lbo etary cia 
1856.14 
AUB oee 
2520, 50 
£322,601 
1609.98 
872.07 
220067 © 
Zo50e2e 


h 


(cm) 


449.13 
908.97 
O34) 
617.22 
666.62 
444.00 
SESS 7 
569 .2u 
620. 75 
CGE oS 
460.43 
DO 
ag Parada 
6227, 06 
671.88 
465.17 
921.54 
D7 De 1 
62 672 SD 
674.01 
483.05 
534.88 
578.34 
625.55 
675.98 
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Table A.1.13 - Dimensions of Member #5 


2, 0)) 


(cm) 


Zee 
249 
21 & 
2aI8 
3.20 
Zee 
2.5L 
ZO 
2599 
SA Alt 
Lie 27 
Zee 
2.16 
3201 
Lae 
PLAS, 
Zeid 


Zia) Oe 


3.02 
Se25 
2.57 
Za 
2.80 
SE 08. 
3.24 


A 
(cm2) 


lies Ls 
194,29 
23967 Ek 
286.00 
553s 29 
134.86 
i972 86 
243.00 
289.00 
336.00 
139.00 
201.14 
246.00 
Pigs SM cal 
338.57 
Owe Ze 
204.00 
249.07 
294.14 
54.0574: 
175.00 
21a. od 
250.86 
296.29 
S426 7 L 


an 





Frame 


(S2 - S5) 


10-12 
10-14 
10-16 
10-18 
10-20 
11-12 
11-14 
11-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
13-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


Table A.1.14 - Dimensions of Member 
2210) 


(Factor Safety = 


Am h 

(cm?) (cm) 

459.76 253.81 
Ww tt 

410.82 239.08 
7 ii] 

364.89 27a a3 

320.90 209.51 
Ww ? 


279293 194.67 
" w 


~9]- 








Frame 


(S2 =- S5) 


10-12 
10-14 
10-16 
LO—iS 
10-20 
eZ 
11-14 
11-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
13-12 
13-14 
3-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


Table A.1.15 - Dimensions of Member 
20) 


(Factor Safety = 


A h 
(cmn2) (cm) 
445.95 249.74 
495.84 264.18 
545, 57 Zi eo 
W wT 
594.13 290.64 


642.58 BUZZ a0 
ul Ww 


oo = 





£22 





Frame 


(S2 = §5) 


10-12 
10-14 
10-16 
10-18 
10-20 
11-12 
11-14 
11-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
L3-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14-20 


(Factor Safety 


A 
(cih2) 


326.590 
265.10 
206. 39 
152.42 
101.89 
326.50 
265510 
206.39 
152.42 
101.89 
326,90 
265.10 
206259 
152.42 
101.89 
326.50 
265.10 
206.39 
152.42 
LOL. 89 
326.50 
265.10 
206.39 
152.42 
LOLs 39 


ao 3— 


Table A.1.16 -—- Dimensions of Member #8 
2.0) 


(cm) 


12 8 
IEE ky; 
1.04 

ook 


1. 28 
diay Ey 
1.04 

Rok 


l. 20 
areal 
1.04 

mil 

aati 
2° 
lla 


sO es 


aod 
Pr i 
Ps: 
a? 
1.04 
oe 
Aa) ae 


A 
F 
(cm? 


) 


he rae fe, 
Aus o 
36.42 
724 ee Ae 
Gis 6S 
Spe a, 
47.79 
36.42 
26 oto. 
LOs63 
9s i 
A779 
36.42 
26.051 
64:63 
Saree, 
47.79 
36.42 
26,0 
ioe o3 
39.50 
re) 
36.42 
ZO el: 
1B) Bs 





Frame 


(S2 = S5) 


10-12 
10-14 
10-16 
10-18 
10=20 
11-12 
11-14 
11-16 
11-18 
11-20 
12-12 
12-14 
12-16 
12-18 
12-20 
13-12 
13-14 
13-16 
13-18 
13-20 
14-12 
14-14 
14-16 
14-18 
14=20 


Table A.1.17 - Dimensions of Member #9 


(Factor Safety = 2.0) 


A 
(om2) 


1529 202 
1240.44 
260529 
698,25 
455.88 
15S 90 2 
1240.44 
960.29 
696525 
455.88 
13S 902 
1240.44 


500.09 
447.21 
391.4 3 
Saad 
264.41 
900.09 
447.21 
391.43 
66) Bae at 
264.41 
200509 
447.21 
391.43 
Silk 27 
264.41 
900.09 
447.21 
391.43 
SS 47) 
264.41 
500.09 
447.21 
391.43 
53L.27 
264.41 
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2.16 
1.86 


2.44 
Zi 
1.86 
1s) 
ZnO 
2.44 


26s, 


1.86 
2 
2.70 
2.44 
Ze LG 
1.86 
2 


A 
(cm2) 


87657 
Tou. 00 
114.91 
82450 
52.43 
Les? 
150.00 
114.91 
82.30 
ara: 
VOt eos 
50-60 
4 
32.250 
52.43 
B75) 
L502 00 
114.91 
82.30 
32 ea 
137 a7 
130,00 
4a 
82.30 
52543 





Table A.1.18 


Limit Design Frame Weights for Single Cross-Tie Tanker 


-95- 


Frame Weight Frame Weight 
Frame (metric tons) (metric tons) 
(S2 - §5) FS = 1.5 PSG= 72 30 
10-12 Ee ee 48.89 
10-14 a2 49.09 
10-16 44.82 oy eae 
10-18 48.96 56.44 
10-20 56.08 63.45 
11-12 41.55 47.85 
11-14 41.99 48.63 
11-16 44.41 oy Gee ee! 
11-18 48.67 55.94 
11-20 oar G2 
12-12 41.47 48.41 
12-14 41.90 48.50 
12-16 44.29 50.96 
12-18 428.53 os ee, 
12-20 on 62.48 
13-12 Ql,72 48.65 
13-14 42.14 48.70 
13-16 44,50 51.09 
13-18 48.68 Bored L 
13-20 5273 62.34 
14-12 41.67 437.51. 
14-14 42.89 48.40 
14-16 44.92 51.54 
14-18 49.11 One 
14-20 5509 62.85 





APPENDIX B - TANKER FRAME WITH RIGID CORNERS 


In large ships which have relatively simple 
web frame geometries, the webs of the individual beams 
are likely to be quite deep. Brackets are also rather 
large. 

In this study, one such frame is examined 
using the assumption of rigid extensions in the bracket 
and beam overlap regions. The web frame numbered 
0C22,49.2t was selected from Lund's paper [5] as a 
representative frame from which a model could be 
developed. This frame was among the lighter frames 
having zero cross<ties and is shown in Figure 10 of 
Reference [5]. 

The frame model is shown in Figure A.3. The 
dots on the system lines represent the end points of 
the rigid extensions. Here, as in Appendix A., the 
system lines are assumed to be defined on the inside 
surface of the shell plating and on the center tank 
face of the longitudinal bulkhead plating. Plastic 
neutral axes of the individual beam elements calculated 
did in fact lie on or very near these lines. The 
principal dimensions of the ship under consideration are 
the same as those listed in Appendix A. 

The loading conditions considered are those 


defined by Lund [5], discussed in Appendix A and shown 


=96= 





in Figure A.2. 

The assumed lengths of bracket legs are 
identical to those used by Lund and are listed in Table 
2 of Reference [5]. 

The lengths of the rigid extensions were 
calculated by hand using the procedures outlined in 
Section 2.2 and Reference [6]. Although Lund only 
calculates these values once for each frame, he does 
not list his results. It is therefore assumed, but 
not known as fact that rigid extension lengths used in this 
study are approximately equal to values used by Lund for 
the same frame. The actual values used are listed in 
Table A.2.1. 

In this analysis the calculations required 
to develop the limit design problem were done by hand. 
The assumption of rigid joint regions required that 
each proposed collapse mechanism be carefully inspected 
in insure kinematic admissibility. Joint rotations 
were not considered as fundamental mechanisms in the 
formations of mechanisms which involved more than one 
frame element. A discussion of the joint rotation 
problem is given in Section 2.3. 

The results of the hand calculations are 
listed in Table A.2.2 The values of the constraint 
Matrix coefficients (i.e. internal energy dissipation 


coefficients) as well as the values of the constraint 
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equation right hand sides (i.e. external work rates) 
are tabulated. 

It should be noted that although the values 
of constraint relations are shown in the table to be 
segregated by load condition, they were combined as one 
system in the linear programming problem. This allowed 
the system to be truly optimized for all considered 
loadings. In general, this is not possible in the 
elastic design methods. Most elastic procedures optimize 
for each load condition and then select the larger 
dimensions required for each element. 

The linear programming subroutine M.P.S.X. 
is used to optimize the results of the hand calculations 
as discussed in Appendix C. Again, asS in Appendix A, 
the output from the procedure is a listing of minimum 
values of plastic limit moment required to make the 
structure safe under the prescribed loading conditions. 

The dimensions of the beam elements were 
calculated based on the relation 

M, = oy Zp 
where 


Z = Aa 
ae 


In these relationships M, is the plastic limit 
moment calculated in the limit design problen, ee is the 


yield stress of the material, A is the total required area 


of the cross section and ais the distance between 
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centroids of the two equal areas lying above and below the 
plastic neutral axis. 
In order to allow comparison to the parent 
Lund frame, the beam elements are sized to meet the 
plastic moment requirements based upon the proportions 
of the beams which Lund used as starting values for 
his frame iterations. For this frame, starting values of 
web depth, web thickness and top flange area are shown 
in Table A.2.3. Also shown are values of plate area 
which are assumed to remain constant. 
Values of beam geometry were calculated using 
FACLOYrS OF Safety Of 1.0, 1.5, 2.0 and 2.5. A corrosion 
allowance of 0.2 cm was applied. Table A.2.4 gives 
resulting values of net beam cross sectional area Ans 
web depth h, web thickness t and top flange area Ap. 
Frame weights were calculated for one half of 
complete frames for each of the factors of safety 


considered. The relation 
6 
= A} = . + 
Ww pt i Ap ~ % (oc, + Ap) 


was used to calculate these weight values where o is 
material density, 1. is the member length and c. and Ch 
are the overlaps at the respective beam ends. This 
expression was also used by Lund. Consistent with Lund's 


calculations, weight of web stiffeners and brackets are 


neglected. 
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Values of frame weight are plotted against 
factors of safety for this analysis and are shown in 
Figure A.6. Also plotted is the weight of the Lund 


elastic design frame at a factor of safety of 1.5. 
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Table A.2.1 - Rigid Extension Lengths 


Rigid Length Rage Length 
Extension (cm) Extension (cm) 
Lp 504 Lu 453 
Le 630 Lay 453 
Ly 128 Lo 153 
Ly 428 Ly 392 
Lp 553 Lo 382 
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Table A.2.3 - Initial Member Dimensions 





Web Top Flange Plate Flange 
Member Height Thickness Area Area 
No. (cm) (cm) (cm) (cm2) 
1 400 2ai0 120 750 
2 400 V8 720 750 
3 300 Pes 120 660 
4 300 Lise 120 420 
5 300 1.4 80 150 
6 240 15, 80 750 
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Table A.2.4 - Dimensions of Limit Design Results 
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APPENDIX C - NUMERICAL ROUTINE 


C.l Use of the M.P.S.X. 

The M.P.S.X. linear programming subroutine used 
in this analysis was found to be very sensitive to format 
errors. The user's guides and instructions available on 
the system are not very careful to point out possible 
pit falls. It is the purpose of this section to list 
the control language statements used in solving the 
linear programming problems of this analysis and to 
point out general problems that were encountered in 
using the M.P.S.X.. Before attempting to use this 
subroutine the user should become very familiar with 
References [31 - 33]. 

A sample listing of the control program with 
sample data points follows: 


//LP EXEC MPSX 

//MPSCOMP.SYSIN DD * 
PROGRAM 
INITIALZ 
MOVE (XDATA, 'LUND') 
MOVE (XPBNAME, 'PBFILE') 
CONVERT (' SUMMARY ' ) 
BCDOUT 
SETUP 
MOVE (XOBJ, 'MINWT' ) 
MOVE (XRHS, 'WORK') 
CRASH 
PRIMAL 
SOLUTION 
RANGE 
EXIT 
PEND 

* 


//MPSEXEC.SYSIN DD * 
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NAME LUND 


MOM1 MECHA 4.0 MECHB 8, 


WORK MECHA 3.605E+8 MECHB 4.228E+7 


ENDATA 
/* 


Explanations of the control language statements 
can be found in Reference [33]. The format sensitivity 
of the M.P.S.X. has been mentioned previously. The 
greatest difficulty occurs when entering data. Page 
209 of Reference [33] will be found to be invaluable 
to any prospective system user. This page tabulates all 
required data entries in the positions required by the 
M.P.S.X.. The data set name and all data vector names 
must be left justified in their fields. Failure to 
have the data set name left justified will result in 
having the data set not read. If a row or vector name 
is not left justified, its associated data pointe will 
not be read. If a row or vector name extends beyond the 
left bound of its field, the data will be read and labeled 
with an erroneous name. Any of these errors associated 
with row and vector names will result in erroneous 
solutions. 

The scale values of the data may be placed 


anywhere in the prescribed field. 
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Finally, Chapter 4 of Reference [33] should be 
checked to insure that data and vector names use only 
"legal" characters. The use of blanks and certain other 
characters will result in abandoned runs. 

Giz Mor. 5<%.- Costs 

The two positive advantages of the M.P.S.X. 
are its power and efficiency. The linear programming 
problems solved in this analysis are moderately large, 
having roughly fifty constraint relations and nine 
variables. There was originally some fear that the 
solutions of these problems would be quite expensive. 

It was found that problems having approximately 
50 constraints and 6 variables required approximately 
0.325 minutes in total job run time at a cost of 
approximately four dollars. These numbers are based 
on computation performed in solving the linear programming 
problem of Appendix B. 

Run time and associated costs were found to be 
very closely related to the total number of constraints. 
Changes in the number of solution variables had little or 
no effect on computation time. 

The total cost of running the limit design 
program in Appendix A for twenty-five different frames 
including the card punching cost was $17.57. The cost 
of solving the programming problems formed by the FORTRAN 


program was $74.08. This resulted in a total computation 
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cost of $91.65 for limit design problem solutions of 


twenty-five frame configurations. 
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